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Abstract X-ray observations have played a key role in the study of substruc- 
ture and merging in galaxy clusters. I review the evidence for cluster 
substructure and mergers obtained from X-ray observations with satel- 
lites that operated before Chandra and XMM. Different techniques to 
study cluster mergers via X-ray imaging and spectral data are discussed 
with an emphasis on the quantitative analysis of cluster morphologies. 
I discuss the implications of measurements of cluster morphologies for 
cosmology and the origin of radio halos. 
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1. INTRODUCTION 

"Substructure" in a galaxy cluster is defined as multiple peaks in the 
cluster surface density on scales larger than the constituent galaxies; the 
"cluster surface density" refers either to the galaxies, the X-ray emission 
from hot gas, or the dark matter. Today we take it for granted that 
many galaxy clusters exhibit substructure and thus are in early stages of 
formation. This, of course, was not always the case. In the 1980s there 
were several searches for cluster substructure in the optical, but their 
results were controversial, primarily because of the difficulty in assessing 
the importance of projection effects and the statistical significance of 
substructure (see reviews by West 1990, 1995). 
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X-ray studies of clusters are less susceptible to contamination from 
foreground and background objects than optical studies. The X-ray lu- 
minosity is a strong function of the temperature, or mass, which means 
that, e.g., foreground groups contribute proportionally less to the X-ray 
emission than they do to the galaxy surface density. X-ray studies of 
clusters also have the advantage that the signal is limited only by the ef- 
fective area of the detector and exposure time of an observation whereas 
optical studies are limited by the finite number of cluster galaxies. 

The reality of substructure in clusters was firmly established with 
ROSAT observations in the early 1990s. The watershed example is that 
of A2256 which had long been thought to be a prototypical relaxed clus- 
ter when examined from the perspective of its galaxy isopleths. However, 
in a controversial optical study of A2256, Fabricant, Kent, and Kurtz 
(1989) proposed the existence of an infalling subcluster from analysis of 
the galaxy velocities even though they detected no such evidence from 
the galaxy positions alone. 




Figure 1.1. (Left) ROSAT PSPC image of A2256 (Briel et al. 1991). (Right) 
ROSAT PSPC image of Coma (Briel & Hemy 1997). 

The existence of a subcluster in A2256 was confirmed by the stunning 
ROSAT PSPC image (Briel et al. 1991) that showed a subcluster^ offset 
from the main cluster by a few hundred kpc (Figure 1.1, left). ROSAT 
also clearly demonstrated significant subclustering in the Coma cluster 
(Briel, Henry, & Bohringer 1992) which had been presumed to be the 
quintessential relaxed cluster (Figure 1.1, right). Hence, ROSAT im- 
ages confirmed and clearly established the existence of substructure in 
clusters, and thus showed that such clusters are really still forming. 
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Figure 1.2. Contour plots of ROSAT PSPC images (see Buote & Tsai 1995) of 
four Abell clusters labeled according to the Jones & Forman (1992) morphological 
classification scheme. 



The fundamental question raised by these early ROSAT observations 
is how widespread is merging in clusters? Are clusters generally young 
or old? Or is there an equal distribution of cluster ages in a given cluster 
sample? To address this issue one needs to have measurements of the 
sub clustering properties of a large cluster sample and, of equal impor- 
tance, a precise definition of the "age" of a cluster. The first systematic 
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X-ray study of cluster merging was by Jones & Forman (1992). From 
visual inspection of ~ 200 Einstein cluster images, Jones & Forman sep- 
arated the clusters into 6 morphological classes (see Figure 1.2). These 
classes range from relaxed single-component systems to systems with 
a large degree of substructure. From the relative populations of these 
classes they deduced that ~ 30% of clusters have substructure, which is 
actually a lower limit because of the limited resolution of the Einstein 
IPC. This study established that merging and substructure are very 
common in clusters. Consequently, the need arose for a more precise 
assignment of the age of a cluster; e.g., how much older or younger are 
clusters in the Jones &; Forman classes? Hence, Jones &; Forman (1992) 
ushered in the era of quantitative X-ray cluster morphology. 

2. QUANTITATIVE ANALYSIS OF 
INDIVIDUAL SUBSTRUCTURES 

Quantitative studies of cluster X-ray morphologies have traveled down 
two distinctly different paths. The first path is that of the detailed 
structural analysis of clusters to determine the number of substructures, 
their fluxes, spatial properties, etc.. A popular approach is to examine 
the residuals obtained from subtracting a smooth model representing a 
relaxed cluster from the X-ray cluster image (e.g., Davis &; Mushotzky 
1993; White et al. 1994; Davis 1994; Prestwich et al. 1995; Neumann & 
Bohringer 1997, 1999; Bohringer et al. 2000). Usually this smooth model 
is obtained by fitting a set of perfect elliptical isophotes or an elliptical 
(3 model to the cluster surface brightness; i.e., the X-ray emission of a 
relaxed cluster is assumed to be elliptical in shape. 

In hydrostatic equilibrium the surfaces of constant X-ray emissivity 
are identical in shape to the surfaces of constant gravitational poten- 
tial regardless of the temperature profile of the gas (Buote & Canizares 
1994, 1998). And since the isopotential surfaces generated by an ellipti- 
cal matter distribution (which is assumed to be the most general stable, 
relaxed, non-rotating, self-gravitating configuration) are not perfect el- 
lipses (e.g., Binney & Tremainc 1987), neither are the X-ray isophotes. 
Consequently, the residuals obtained from subtracting elliptical mod- 
els from the X-ray surface brightness of clusters need to be carefully 
considered. This procedure is most appropriately applied as a simple, 
approximate indicator of substructure. 

A more general and powerful method to identify and quantify sub- 
structures is to perform a wavelet decomposition of the X-ray image. 
The wavelet analysis is a powerful multi-scale technique to detect sources 
embedded in the bright diffuse background cluster emission which has 
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been successfully applied to many clusters (e.g., Slezak et al. 1994; 
Vikhlinin, Forman, & Jones 1994; Grebenev et al. 1995; Biviano et al. 
1996; Pislar et al. 1997; Lima-Neto et al. 1997; Pierre & Starck 1998; 
Lemonon et al. 1997; Dantas et al. 1997; Vrtlik et al. 1997; Lazzati & 
Chincarini 1998; Lazzati et al. 1998; Arnaud et al. 2000). Wavelet anal- 
ysis locates substructures on different scales and allows separate spatial 
analysis (e.g., flux, extent etc.) of each detected structure. The sta- 
tistical significance of the substructures can be assessed rigorously via 
Monte Carlo simulations. 




Figure 1.3. Wavelet decomposition of ROSAT images of (Left) A2256 by Slezak et 
al. (1994) and (Right) of Coma by Biviano et al. (1996). 



Applications of wavelets to the ROSAT images of A2256 and Coma 
(Figure 1.1 are shown in Figure 1.3. In the case of A2256 Slezak et 
al. (1994) establish that the core is more than a simple bimodal system 
since the bottom-left region consists of at three subclusters The wavelet 
analysis of Coma by Biviano et al. (1996) shows that the core consists 
of two subclusters surrounding each of the large galaxies NGC 4874 and 
NGC 4889. Apparently both Coma and A2256 are far from relaxed 
systems. 

Wavelets are particularly useful for less-massive systems like A1367 
where the emission from several galaxies or groups needs to be separated 
from the diffuse cluster background. In Figure 1.4 is shown the wavelet 
analysis by Grebenev et al. (1995) who analyzed both the ROSAT PSPC 
and HRI images and detected 16 extended sources embedded in the 
diffuse ICM of A1367. Not only does the wavelet analysis allow the 
fluxes and extents of each of these sources to be quantified, but the 
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Figure 1.4- ROSAT image of A1367 and a wavelet decomposition on different scales 
by Grebenev et al. (1995). 



larger scale wavelets (see Figure 1.4) show that the cluster is bimodal 
with subclusters centered about what are likely to be galaxy groups. 

The power of the wavelet technique is also demonstrated by the anal- 
ysis of the ROSAT HRI image of A521 by Arnaud et al. (2000). From 
visual inspection of the HRI image of A521 one notices asymmetric 
isophotes such that the emission peak appears to be offset from the 
centers of the fainter isophotes (see left panel in Figure 1.5). Applica- 
tion of the wavelet technique to this image reveals two distinct structures 
(see right panel of Figure 1.5). The main cluster appears to be oriented 
along a line connecting two adjacent clusters. Nearly perpendicular to 
this line is the line connecting the subcluster to the main cluster. This 
other line appears to lie nearly parallel to the line pointing to another 
adjacent cluster. Consequently, Arnaud et al. conjecture that A521 lies 
at the intersection of two large-scale filaments. 
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Figure 1.5. (Left) ROSAT HRI image of A521 and (Right) the wavelet transformed 
image (Arnaud et al. 2000). 



3. QUANTITATIVE CLASSIFICATION OF 
GLOBAL MORPHOLOGY 

The other path taken by studies of quantitative X-ray cluster mor- 
phology is to build on the work of Jones &: Forman (1992) and to devise 
a quantitative scheme for classifying the morphologies of X-ray images 
of galaxy clusters. As with any classification system in astronomy the 
principal motivation for classifying cluster morphologies is to elucidate 
fundamental physical properties, in particular those associated with clus- 
ter formation and evolution. 

The presence of substructure in clusters implies they are still forming 
and evolving dynamically, and thus a logical candidate for a fundamen- 
tal parameter is the current dynamical state. The dynamical state of 
a cluster is related to the amount of time required for the cluster to 
virialize; i.e., a time of order a crossing time. But for a cluster of a 
given total mass one can imagine many different morphological configu- 
rations - and formation histories - that would lead to similar relaxation 
timescales. Hence, to classify clusters having different formation histo- 
ries but similar dynamical states we also require one or more fundamen- 
tal parameters to specify the type of merger (e.g., bimodal, many small 
subclusters) as indicated qualitatively by the classes of Jones &: Forman 
(1992). 
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3.1. METHODS 

Perhaps the most common approach used to quantify the morpholo- 
gies of a large number of X-ray cluster images has been with a measure 
of the X-ray elhpticity (e.g., McMillan et al. 1989; Davis 1995; Mohr 
et al. 1995; Gomez et al. 1997; Gomez, Hughes, & Birkinshaw 2001; 
Kolokotronis et al. 2001). This method is not a particularly good indi- 
cator of the dynamical state since both relaxed and disturbed clusters 
can have significant cllipticity. And even disturbed chistcrs can have 
small ellipticity if the substructure is distributed symmetrically about 
the cluster center. Moreover, even if both the ellipticity and associated 
position angles are considered they only provide a crude measurement of 
cluster morphology and have never been shown to provide an interesting 
distinction between the variety of morphologies exemplified by the Jones 
k, For man classes. 

A better method is the center-shift introduced by Mohr, Fabricant, & 
Geller (1993). This popular method has been applied in various forms 
to X-ray cluster images in several studies (e.g., Mohr et al. 1995; Gomez 
et al. 1997, 2000; Rizza et al. 1998; Kolokotronis et al. 2001). The 
basic idea is to divide up a cluster image into a series of circular annuli 
having different radii but with centers located initially at a guess for the 
cluster center. The center-shift is then given by a weighted average of 
the centroid computed for each of these annuli. 

Since the center-shift is sensitive only to asymmetries in the X-ray 
images (in particular non-ellipsoidal configurations) it is a much more 
reliable than the cllipticity as an indicator for when a cluster is relaxed. 
However, it is not transparent how the center shift translates into a 
physical measure of the dynamical state. And since the center-shift is 
most sensitive to mergers of equal-mass subclusters, it cannot by itself 
distinguish the full range of structures exhibited by the Jones & Forman 
morphological classes. 

If the only objective were to distinguish the full range of cluster mor- 
phologies then the logical procedure would be to decompose cluster im- 
ages into a set of orthogonal basis functions of which wavelets (see §2) are 
the probably best example. The wavelet coefficients would then define 
the parameter space of cluster morphologies. Unfortunately, there is no 
obvious connection (of which I am aware) between wavelet coefficients 
and a physical measure of the dynamical state. 

One method that is both closely related to the cluster dynamical state 
and provides a quantitative description of the full range of Jones & 
Forman morphological classes is the "power ratio" method (Buote & 
Tsai 1995, 1996; Buote 1998). The power ratios are constructed from 
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the moments of the two-dimensional gravitational potential. Specifically, 
one evaluates the square of the moments over a circle of radius, R, where 
the origin is located at the center of mass or the at the largest mass peak. 
The ratio of term, m, to the monopole term is called a "power ratio" , 

Po - [in'?) ' ^ ^ 

where "f™* is the mth multipole of the two-dimensional gravitational 
potential due to matter interior to the circle of radius, i?, and (• • •) 
represents the azimuthal average around the circle. In detail we have, 

Po = [aoln(i?)]', (2) 

for m = 0, 

' -.2 , 1-2 



a„ + 6„ (3) 



for m > 0. The moments am and hm are given by, 

am{R) = I Y.{x') iR')"^ cos met)' cPx', 

Jr'<r 

bm{R) = [ J:{x')(R')"'smm(p'd^x', 
Jr'<r 

where x' = {R', (f)'). 

These ratios are directly related to the 2D gravitational potential if 
one has a map of the 2D surface mass density such as provided by weak 
gravitational lensing studies. For X-ray studies S is replaced with the 
X-ray surface brightness. Ex, and therefore the power ratios in X-ray 
studies arc really derived from a pscudo potential. These ratios are 
most sensitive to structures on the same scale as the aperture radius, R. 

When the aperture is located at the peak of the X-ray emission the 
dipole power ratio, Pi/Pq, provides structural information similar to 
the center shift discussed above (see also Dutta 1995). For an aperture 
located at the centroid of the surface brightness the dipole moment van- 
ishes. In this case the quadrupole power ratio, P2/P0, is sensitive to the 
degree of flattening and is related to the ellipticity. But unlike ellipticity 
P2/P0 is also sensitive to the radial profile of the X-ray emission. 

The primary physical motivation behind the power ratios is that they 
are related to potential fluctuations. And since it is thought that large 
potential fluctuations drive violent relaxation in clusters, the power ra- 
tios are closely related to the dynamical state of a cluster (Buote 1998). 
The other motivation is that the multipoles are a complete orthogonal 
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set of basis functions for the (pseudo) potential and thus are well suited 
to classify the wide range of observed cluster morphologies. 

To get a feel for the power ratios let us see how they behave on the 
ROSAT PSPC images of clusters in the different Jones & Forman mor- 
phological classes shown in Figure 1.2. The four clusters inhabit the 
extreme Jones & Forman classes. A2029 is a smooth, single component 
system apparently in a relaxed state. A85 has a regular dominant com- 
ponent but with a small structure ~ 0.6 Mpc to the S. A1750 is a double 
cluster consisting of two roughly equal-sized components separated by 
~ 1 Mpc. A514 is a highly irregular aggregation of structures. 




Figure 1.6. Power ratios (from Buote & Tsai 1996) for the clusters in Figure 1.2 
computed within a circular aperture of 1 Mpc radius located at the centroid of the 
X-ray emission. 



In Figure 1.6 I show the power ratios, P2/P0 and P3/P0, of these 
clusters computed for a 1 Mpc aperture^ where the aperture is located 
at the centroid of the X-ray emission (i.e., analog of the center of mass). 
It can be seen that the single-component cluster is well separated from 
the primary with small secondary. And each of these classes is clearly 
distinguished from the disturbed complex and double clusters. In effect 
the power ratios have defined a morphological evolutionary track where 
the young, unrelaxed clusters are born at the top right of the figure. 
As they relax and erase their substructure they pass through a phase 
similar to A85 until they are old and evolved systems like A2029. 
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Although we have succeeded in obtaining a successful broad classifi- 
cation according to dynamical states, we have not distinguished clearly 
between the different classes of highly disturbed clusters (i.e., complex 
and double). Since there is nothing special about the 1 Mpc aperture it 
is sensible to explore the effects of using different apertures. The result 
of computing the power ratios in a 0.5 Mpc aperture are displayed in 
Figure 1.7 
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Figure 1.7. As Figure 1.6 but for the 0.5 Mpc aperture. 



By focusing initially on P2/P0 it can be seen that three of the clusters 
appear to be relaxed systems (i.e., small P2/Po)- This is because the 0.5 
Mpc aperture only encloses 1 component of the double cluster and only 
the primary component of A85. The single component cluster A2029 
appears relaxed on both the 0.5 and 1 Mpc scales. However, A514 is 
complex on many scales and it is easily distinguished from the other 
reference clusters as a disturbed system in the 0.5 Mpc aperture. Of 
course, one only needs to appeal to P3/P0 to verify that both the double 
and complex clusters are actually in a younger dynamical state than 
the others. Hence, the power ratios represent a quantitative 
implementation of the Jones & Forman classification scheme, 
particularly on the 0.5 Mpc scale. 
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Figure 1.8. Power ratios of the brightest ~ 40 clusters (Buote & Tsai 1996) com- 
puted within apertures of 0.5 Mpc (Left) and 1 Mpc (Right). 



The result of computing power ratios for the brightest ~ 40 ROSAT 
clusters is displayed in Figure 1.8. It is immediately apparent that there 
is a marked deficiency of highly disturbed clusters (complex and double). 
These brightest clusters therefore lack young members and are instead 
dominated by mostly evolved clusters with only small-scale (< 500 kpc) 
substructure. Since such highly evolved clusters are usually associated 
with cooling flows it should be expected that cooling flows dominate 
the brightest clusters as has been suggested before on different grounds 
(e.g., Arnaud 1988; Forman & Jones 1990; Edge et al. 1992; Peres et al. 
1998). 

In Figure 1.9 the quantitative connection between cooling flows and 
cluster morphology is shown by the anti-correlation of the mass deposi- 
tion rate (M) and P2/P0. This represents the first quantitative descrip- 
tion of the anti-correlation of substructure with the strength of a cooling 
flow. Note the large scatter for systems that have significant substruc- 
ture (i.e., large P2/ Pq). Analysis of this correlation and its large scatter 
should shed light on how cooling flows are disrupted by merges and are 
subsequently re-established. 
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Figure 1.9. As Figure 1.8 for the 1 Mpc aperture except that the coohng flow mass 
deposition rate has been plotted on the vertical axis. 
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Figure 1.10. ROSAT HRI images of (Left) RXJ 1347- 1145 from Schindler et al. 
(1997) and (Right) CL0024+17 from Bohringer et al. (1997). 



4. HIGH-REDSHIFT CLUSTERS 

Unfortunately, because of the limited resolution and collecting area 
of ROSAT it has been difficult to study the morphologies of distant 
clusters. Two of the best examples {z ~ 0.4) imaged with the ROSAT 
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HRI are displayed in Figure 1.10. The cluster RXJ1347.5-1145 appears 
to be a relaxed, cooling flow (Schindler et al. 1997) while the cluster 
C10024+17 may have substantial substructure as quantified by a center 
shift (Bohringer et al. 2000). These tantalizing glimpses demonstrate 
the need for a systematic study at high resolution with Chandra. 

5. MORPHOLOGY AND COSMOLOGY 

Fossil imprints of the process of formation arc retained in the cluster 
substructure. In the standard hierarchical paradigm of structure forma- 
tion the mass spectrum of subclusters is related to the power spectrum of 
mass density fluctuations which is a key distinguishing property of cos- 
mological models (e.g., Peacock 1999). As clusters evolve dynamically 
the mass spectrum of subclusters changes. In a standard Friedmann- 
Robertson- Walker universe with < 1 and Aq = 0, the linear growth 
of density fluctuations becomes strongly suppressed when the curvature 
term in the Friedmann equation exceeds the matter term. The redshift 
delineating this transition from an Einstein - de Sitter phase to one of 
free expansion is then 1 + Ztrans = — 1; i-e. when the matter density 
f2(ztrans) = 0.5. Hencc, if r^o ^ li then objects formed a long time in the 
past relative to universes where ~ 1, and thus clusters in low-density 
universes should be, on average, more relaxed than clusters in universes 
with J7o ~ 1. 

5.1. SEMI-ANALYTICAL MODELS 

Richstone, Loeb, & Turner (1992) presented the first theoretical model 
relating i7o to the observed frequency of substructure in clusters. In 
their semi-analytical calculations they avoided the issue of the power 
spectrum by concentrating on clusters having the same total mass. The 
collapse time of a 10^^ /i~^Mq spherical density perturbation (taken to 
be twice the turn around time) was defined to be the dividing point 
between clusters that do and do not possess substructure. By further 
assuming that any substructure is erased on a crossing time (taken to 
be 0.1/Hq), Richstone et al. computed the quantity 5F, the fraction of 
present-day clusters which formed within the last time interval, 0.1 /Hq, 
as a function of and Aq. They found dF ~ Qq (see Figure 1.11). When 
compared to the estimates of ^ 30% for the frequency of substructure 
in nearby clusters (Jones &; Forman 1992) Richstone et al. concluded 
that J^o ^ 0.5. 

Follow up theoretical studies by KauflFmann Sz White (1993), Lacey 
Sz Cole (1993), and Nakamura, Hattori, & Mineshige (1995) emphasized 
that the time for substructure to be erased is variable and can be es- 
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Figure 1.11. Figure 3 from Richstone et al. (1992). Fraction of present-day clusters 
which formed within the last time interval, St, as a function of flo. The time inter- 
vals are in units of Hq^. The horizontal line represents a conservative observational 
estimate of the frequency of substructure. 

pecially long for substructures with compact cores. The relationship 
between the collapse time of a spherical density perturbation and sub- 
clustering, though qualitatively reasonable, is ambiguous. Consequently, 
it is difficult to compare directly the frequency of observed substructure 
to predictions of semi-analytic models based on Richstone et al.'s idea. 

Thus, a fundamental limitation of these studies is that they only 
predict the ambiguous "frequency of substructure" rather than a well- 
defined quantitative measure of cluster morphology such as the power 
ratios. Since Richstone et al.'s idea is really a statement about the dy- 
namical states of clusters, in Buote (1998) I used a related (but more 
detailed) semi-analytical approach to study the behavior of cluster power 
ratios in different cosmologies. Violent relaxation (Lynden-Bell 1967) is 
the key process driving the elimination of large potential fluctuations. It 
operates on a timescale of ~ 1 — 2 crossing times and proceeds indepen- 
dently of the masses of the constituents. Consequently, I argued that a 
plausible definition of the dynamical state of a cluster is. 
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where M is the average mass and AM is the mass accreted over a relax- 
ation time, typically assumed to be a crossing time; AM/M is called the 
"fractional accreted mass". This equation states that over the duration 
of a crossing time the fractional increase in the rms spherically averaged 
potential is approximately equal to the fractional increase in the mass 
added in quadrature to the ratios of the increases of the rms spherically 
averaged higher order potential multipolcs to the monopole. 

The key premise is that the amount of accreted mass over the previous 
relaxation timescale determines the amount of substructure (or non- 
ellipsoidal distortions) which is similar to the premise of Richstone et 
al. that substructure is related to the collapse and crossing times. This 
premise requires that AM/ M be strongly correlated with the other low- 
order terms, which are approximately, ((^*J"'*)^)/((^*cr*)^)' defined at 
the epoch of interest. These terms are just the 3D versions of the power 
ratios (see equation 1). 
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Figure 1.12. Prom Buotc (1998): (Left) The mass- averaged fractional accreted mass 
evaluated for r — Ihr^ Mpc at z = for CDM models as a funetion of f2o (Ao = 0). 
The crosses indicate a mass average over the full range (0.35 — 3) x lO^^ft^^MQ, and 
the boxes indicate a lower limit of 7 x lQ^'^h~^ Mq. Relaxation timescales of 1 and 
2 crossing times arc shown. (Right) ^Mj M for clusters of mass 7 x 10^*/i^"'A/q 
evaluated for r = l/i^^ Mpc at « = for models with fio = 1 and power spectra 
P(fc) oc A;" as a function of spectral index n. 

The dependence of AM/ M on fio and the power spectrum is shown 
in Figure 1.12. We see the expected increase in fractional accreted mass 
with increasing f^o where AM/ M oc -v/Hq, but the normalization does 
depend sensitively on the assumed relaxation timescale similarly to the 
previous related studies by Richstone et al. and others. Also shown is the 
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dependence on n, the spectral index of models with P{k) oc k'"', which is 
considerably steeper, AM/ M oc (— ra)^'^. Since the observable low-order 
power ratios should behave as P„JPq ~ (AM/M)^ (see section 2.2 of 
Buote 1998), the power-ratio distribution for a large sample of clusters 
should be an interesting probe of J7o and the power spectrum. 

5.2. N-BODY SIMULATIONS 

N-body simulations of CDM clusters confirm that the mean value of 
Pm/Po for small m in a cluster sample increases with Q,q (Buote & Xu 
1997; Thomas et al. 1998). But Buote & Xu (1997) also perform simu- 
lations with P{k) oc A;"" for different n and find that the mean vahic of 
Pm/Po is barely affected by n. On the other hand they find that n does 
affect significantly the variance of P^/Pq- (^o does not seem to affect 
the variance.) These conclusions have to be viewed with some caution 
because these dark-matter-only simulations analyze the projected square 
of the mass density in an attempt to mimic X-ray observations. Further 
work with large high-resolution N-body simulations is required to estab- 
lish precisely the relationships between AM/ M, P^n/Po, ^o, and P{k) 
(and Ao). 

Other N-body simulations with and without gas show that center-of- 
mass shifts are also sensitive to (Jing et al. 1995; Crone, Evrard, &: 
Richstone 1996). Generally both semi-analytic models and dark-matter- 
only N-body simulations agree that center shifts and power ratios can 
distinguish between CDM models with different values of JIq- The same 
holds for gas-dynamical N-body simulations (Evrard et al. 1993; Mohr 
et al. 1995) 

However, when the N-body simulations (with or without gas) are com- 
pared to X-ray observations of clusters conflicting results are obtained 
(Figure 1.13). Mohr et al. (1995) compare center shifts of clusters 
formed in hydrodynamical simulations to Einstein clusters and conclude 
that J7o ~ 1 whereas Buote Sz Xu (1997) compare power ratios of the 
projected square of the dark matter density to ROSAT clusters and con- 
clude r^o < 1- Furthermore, the clusters formed in the hydrodynamical 
simulations by Valdarnini, Ghizzardi, & Bonometto (1999) give power 
ratios different from those obtained by Buote & Xu (1997). 

All of these simulations have deficiencies. The most important de- 
ficiency in the hydrodynamical simulations is the poor force resolution 
for the gas: softening lengths of ~ 80 kpc for Valdarnini et al. (1999) 
and over 100h~^ kpc for Mohr et al. (1995) . The simulations of Mohr 
et al. also contained only six clusters which is too small for statistical 



18 



• — I 



OCDM (n„ = 0.2. Ao = 0) 



LCDM (fij = 0.2, \ = 0.8) 



SCDM (fio = 1.0, Xo = 0) 
• 1 



EINSTEIN Clusters 



I I 



b 0.6 - 



— 1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I— 



LCDM 



SCDM 



OCDM 



ROSAT 



BCDM 



I ■ 



20 40 60 80 

Center Shift (hgjkpc) 



100 



—I I I I I I I I I I I I I I I l_ 



-5.8 -6.6 -5.4 -5.2 

Average PRj 



Figure 1,13. (Left) Center shifts obtained by Mohr ct al. (1995) from simulated 
and Einstein clusters. (Right) Power ratios {PR2 = logio(-P2/^'o)) obtained by Buote 
& Xu (1997) for ROSAT and simulated clusters: OCDM (Ho = 0.35), LCDM (Qo = 
0.35, Ao = 0.65), and SCDM (^o = 1.0). 



studies. Finally, the simulations of Buote & Xu (1997) approximated 
the gas distribution using the dark matter. 

Clearly until appropriate simulations are applied to this problem we 
will not have a reliable constraint on or P{k) from cluster morpholo- 
gies. What is needed are high-resolution 20 kpc) three-dimensional 
gas-dynamical simulations of a large number 50) of clusters. The 
existing observational samples of Einstein data (Mohr et al. 1995) and 
ROSAT data (Buote & Tsai 1996) also need to be expanded and re- 
analyzed with new high-resolution, high S/N Chandra and XMM data. 
These requirements are not excessive for a problem that deserves serious 
attention. 

6. MORPHOLOGY AND RADIO HALOS 

It has been noticed for some time that X-ray observations provide 
circumstantial evidence for a connection between cluster merging and 
radio halos (see Feretti 2000 and references therein) because, in particu- 
lar, radio halos are only found in clusters possessing X-ray substructure 
and weak (or non-existent) cooling flows. However, it has been argued 
(e.g., Giovannini & Feretti 2000; Liang et al. 2000; Feretti 2000) that 
merging cannot be solely responsible for the formation of radio halos 
because at least 50% of clusters show evidence for X-ray substructure 
(Jones Sz Forman 1999) whereas only ~ 10% possess radio halos. (Note 
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X-ray and optical substructures are well-correlated - Kolokotronis et al. 
2001.) 

Unfortunately, it is difficult to interpret the importance of merging 
using the observed frequency of substructure as it does not itself quan- 
tify the deviation of an individual cluster from a virialized state. And 
the shocks that could be responsible for particle acceleration will be 
proportionally stronger in clusters (of the same mass) with the largest 
departures from a virialized state. To measure the dynamical states of 
clusters from X-ray images it is necessary to quantify the cluster mor- 
phologies using statistics such as the center-shift and the power ratios. 
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Figure 1.14- From Buote (2001): Radio power (Pi. 4 - 1.4 GHz rest frame) versus 
dipole power ratio (Pi/Po) where Pi. 4 includes emission from (Left) only radio halos 
and (Right) the total diffuse emission from halos and relic sources. The power ratios 
are computed within a 0.5 Mpc aperture centered on the X-ray emission peak with 
estimated la errors shown. (Uncertainties on Pi. 4 are believed to be ^ 10% and are 
not shown.) 



In Buote (2001) I used power ratios to provide the first quantitative 
comparison of the dynamical states of clusters possessing radio halos. A 
correlation between the 1.4 GHz power (Pi. 4) of the radio halo (or relic) 
and the magnitude of the dipole power ratio (Pi/Pq) was discovered such 
that approximately Pi. 4 oc Pi/Pq (sec Figure 1.14). The Pi. 4 — Pi/Pq 
correlation not only confirms previous circumstantial evidence relating 
the presence of radio halos to mergers but, more importantly, establishes 
for the first time a quantitative relationship between the "strength" of 
radio halos and relics (Pi. 4) and the "strength" of mergers (Pi/Pq); 
i.e., the strongest radio halos appear only in those clusters currently 
experiencing the largest departures from a virialized state. Moreover, 
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in the Pi a — Pi/Pq plane both radio halos and rehcs may be described 
consistently which provides new evidence that both halos and relics are 
formed via mergers. The Pi. 4 — Pi/Pq correlation supports the idea that 
shocks in the X-ray gas generated by mergers of subclusters accelerate (or 
re-accelerate) the relativistic particles responsible for the radio emission. 

From additional consideration of a small number of highly disturbed 
clusters without radio halos detected at 1.4 GHz, and recalling that radio 
halos are more common in clusters with high X-ray luminosity (Giovan- 
nini, Tordi, & Feretti 1999), I argued that radio halos form preferentially 
in massive {L^ ^ 0.5 x 10^^ erg s~^) clusters experiencing violent merg- 
ers {Pi/Pq ^ 0.5 X 10~^) that have seriously disrupted the cluster core. 
The association of radio halos with massive, large-Pi/Poi core-disrupted 
clusters is able to account for both the vital role of mergers in acceler- 
ating the relativistic particles responsible for the radio emission as well 
as the rare occurrence of radio halos in cluster samples. 

On average Pi/Pq is expected to increase with increasing redshift 
owing to the higher incidence of merging (Buote 1998) which would 
lead to a higher incidence of radio halos. However, on average cluster 
masses are lower at earlier times implying a lower incidence of radio 
halos. Each of these factors is dependent on the assumed cosmology, 
and future theoretical work is therefore required to establish whether 
the abundance of radio halos (1) increases or decreases with redshift, 
and (2) provides an interesting test of cosmological models. 

7. TEMPERATURE SUBSTRUCTURE 

The morphologies of X-ray images of clusters suggest that clusters 
span a wide range of dynamical states and merger configurations. During 
such violent mergers the gas should be shock-heated at various locations 
between an infalling subcluster and the center of the primary cluster. In 
contrast to the azimuthally symmetric temperature profile expected of a 
relaxed system, two-dimensional temperature variations both represent 
a necessary confirmation of the merger picture obtained from images and 
also provide a complementary view of the cluster dynamical state and 
merger history. 

7.1. X-RAY TEMPERATURE MAPS 

In the era before Chandra and XMM it was exceedingly difficult to 
obtain accurate two-dimensional X-ray temperature maps of clusters. 
The ROSAT PSPC had sufficient spatial and spectral resolution but its 
bandpass cut off sharply just beyond 2 keV. Since massive clusters have 
temperatures above ~ 5 keV the temperatures could not be constrained 
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Figure 1.15. ROSAT temperature map of the Coma cluster (Briel & Hemy 1997). 



with any precision for all but a small number of the brightest clusters. 
For these clusters the S/N was so high that the data from the spectra 
below 2 keV managed to place interesting constraints on the tempera- 
ture. 

For example, the ROSAT temperature map of Coma (Briel & Henry 
1997) displayed in Figure 1.15 shows significant temperature variations. 
The region of hotter gas in between the main cluster and the NGC 4839 
subcluster is consistent with shock heating during the passage of the 
subcluster through the main cluster (e.g., Burns et al. 1994; Ishizaka 
&: Mineshige 1996). However, further simulations are required to estab- 
lish whether the subcluster is currently falling in or has already passed 
through the main cluster. As noted by Briel & Henry (1997) if the 
subcluster already passed though the main body then then it is unclear 
why the subcluster still has retained its halo of hot gas. Other ROSAT 
temperature maps of mergers display similar evidence for shock-heating 
(e.g., Briel k Henry 1994; Henry k Briel 1995, 1996; Ettori et al. 2000). 



22 





1 


1 

t 


1 1 1 


1 1 1 1 I 




1 


( ' ' 




Ml,! 

1 1 1 1 





III 


1 


III 

1 


1 

"\\ 


- + 

" 1 


^ \ 

III 


\ 
1 


* t , 

III 


1 



1 2 3 4 5 



1 2 3 4 5 6 7 



10 11 12 13 14 



Figure 1.16. ASCA temperature maps of (Left) Cygnus-A and (Right) A3667 from 
Markevitch et al. (1999). 



The qualitative features in the temperature maps derived from ROSAT 
were also found with data from the ASCA satellite. The higher energy 
resolution and larger bandpass (up to 10 keV) of ASCA provided a dis- 
tinct advantage over ROSAT studies, but the poor spatial resolution 
1.5' FWHM) and highly energy dependent point spread function 
(PSF) seriously hampered two-dimensional spatial-spectral analysis. To 
obtain physical results with ASCA data the PSF needs to be incorpo- 
rated into the analysis. 

When incorporating the PSF into modeling of the ASCA data of merg- 
ers two-dimensional temperature variations similar to those obtained by 
ROSAT are found. For example, in Figure 1.16 the results of the analy- 
sis of Cygnus-A and A3667 by Markevitch, Sarazin, & Vikhlinin (1999) 
are shown. Although some of the detailed results for a particular cluster 
differ between studies using different deconvolution procedures, the basic 
idea that non-azimuthal temperature variations exist in mergers seems 
to be supported by most ASCA and BeppoSAX studies (e.g., Markevitch 
et al. 1998, 1999; Churazov et al. 1999; Donnelly et al. 1999; Molendi 



X-Ray Observations of Cluster Mergers 



23 



et al. 1999; Shibata et al. 1999; de Grandi & Molendi 1999; Henriksen, 
Donnelly, & Davis 2000; Iwasawa et al. 2000). 

Since there are some differences in the radial temperature profiles 
obtained from ASCSA data depending on the PSF deconvolution pro- 
cedure used (see White 1999 and Irwin & Bregman 2000 and references 
therein) the detailed temperature features obtained with ASCA do need 
to be confirmed with Chandra and XMM. (As do those with BeppoSAX 
because of its low spatial resolution.) Nevertheless, the overall trend 
of non-azimuthal temperature structures and the shock-heating of the 
intra-cluster medium are supported by the available ROSAT, ASCA, 
and BeppoSAX data. 

7.2. QUANTITATIVE CLASSIFICATION OF 
TEMPERATURE MORPHOLOGY 
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Figure 1.17. A possible description of the evolution of the X-ray temperature struc- 
ture and image morphology during the formation and evolution of a cluster. 



To obtain a more complete picture of the current dynamical states and 
the merger histories of clusters the global morphological classification of 
cluster images discussed in (§3) should also incorporate the morphologies 
of X-ray temperature maps. In Figure 1.17 I show an idealized picture 
of how the temperature morphology of a cluster might evolve during 
a merger. At early times there is a large spread of temperatures dis- 
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tributed spatially in a non-azimuthally symmetric fashion. At this time 
the cluster is far from a virialized state. It possesses obvious substruc- 
ture and a disturbed spatial morphology quantified by, e.g., a large value 
for the P2/P0 power ratio. The disturbed morphology implies there is 
no cooling flow at this early time (see Figure 1.9), and the cooling time 
(^cooi) is longer than the cluster age (tciustcr) (e-g-, Fabian 1994) 

As the system relaxes image substructure and the spatial fluctuations 
in the temperature are gradually erased until the system is approxi- 
mately isothermal and on the verge of establishing a cooling flow. At 
this time tcool ~ ^cluster a-nd there is only a small amount of substructure 
(i.e., intermediate P2/P0 values - see §3). If the cluster now experi- 
ences a major merger it will begin again at the top of Figure 1.17 with 
a lot of temperature variations and image substructure. If instead the 
system relaxes further without being disturbed then a cooling flow will 
develop (icooi < ^cluster) and the image substructure should be mostly 
erased (smallest values of P2/Po)- Although the azimuthal tempera- 
ture variations will also be erased, a radial temperature gradient will 
be established where the temperature rises from the center out to an 
approximately isothermal plateau. 

Such radial temperature gradients arc characteristic of cooling flows 
(e.g.. White 1999 and references therein). Whether the temperature 
profile is caused by cooling gas or a two-phase medium (e.g., Ikebe et al. 
1997; Xu et al. 1998) is not important for the arguments presented here. 
All that is required is that relaxed systems (particularly those with cD 
galaxies) tend to have characteristic radial temperature structure. 

Therefore, for the merger scenario displayed by Figure 1.17 the amount 
of image substructure (P2/-Fb) fells continuously as the cluster relaxes, 
but the overall spread in temperatures falls and then rises again at late 
times. One possible way to quantify the temperature morphology is 
with the multiphase strength (Buote, Canizares, & Fabian 1999) which 
essentially measures the width of the differential emission measure, ^ (T) , 



where ^max is the maximum value of ^(T) and (T) is the emission- 
measure weighted value of T. This statistic ignores the spatial infor- 
mation and is therefore intended as a relatively crude measure of the 
temperature variations in a cluster appropriate when the data do not 
allow precise temperature estimates in small spatial regions. In such 
cases where the integrated cluster spectrum is modeled with a simple 
cooling flow spectral model plus an isothermal component then equa- 
tion 5 is modified to /a^, where / is the relative fraction of the cooling 





X-Ray Observations of Cluster Mergers 



25 



flow to the total emission measure (see section 5.2 of Buote et al. 1999). 
A variation on this prescription using the breaks in cooling flow mass 
deposition profiles has been used to determine the "ages" of some bright 
coohng flow clusters with ROSAT (Allen et al. 2000). 

Joint consideration of cr^ and Pm/Po should provide a more precise 
indicator of the current cluster dynamical state and merger configuration 
than Pm/Po alone. For high precision temperature maps adding a first 
or second radial moment to equation (5) may be sufficient to capture the 
spatial dependences accounted for in the scenario represented in Figure 
1.17. 

Finally, the scenario described by Figure 1.17 will be complicated if 
there are important dynamical contributions from non-thermal processes 
such as AGN feedback (e.g., Owen ct al. 2001). Empirical studies of the 
spatial and spectral morphologies of a large number of clusters using the 
improved instruments on Chandra and XMM will help to elucidate the 
importance of these and other process associated with cluster formation 
and evolution. 

8. CONCLUSIONS 

X-ray images of clusters obtained by Einstein and ROSAT have es- 
tablished that substructure and merging are common in nearby galaxy 
clusters. This evidence is reinforced by the X-ray temperature maps of a 
smaller number of bright clusters analyzed by ROSAT, ASCA, and Bep- 
poSAX. The study of substructure and morphology has evolved beyond 
detection and visual classification to that of quantitative morphological 
statistics that probe the dynamical states and and the power spectrum 
of density fiuctuations. 

Unfortunately, the present status of cosmological studies of cluster 
morphologies is ambiguous. Although theoretical studies agree that clus- 
ter morphologies are sensitive to the cosmology (particularly to Qq and 
P{k)), the nature of the agreement and the relationship to observations 
have been often in conflict. It is difficult to interpret these disagreements 
because all of the N-body simulations applied to this problem have been 
inadequate. Large volume, high-resolution gas dynamical N-body simu- 
lations are required to obtain definitive answers. A larger observational 
sample of cluster morphologies with higher S/N data is also needed. 

A quantitative connection between cluster mergers and the formation 
of radio halos has now been established. The strength of a merger indi- 
cated by the dipole power ratio (Pi/Pq) is approximately proportional 
to the power of the radio halo. Radio halos form preferentially in merg- 
ers of massive clusters with large values of Pi / Pq where the merger has 
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proceeded fully into the core of the cluster. Larger samples are needed 
to understand the relative importance of the mass and Pi/Pq on the 
strength of the radio halo and to clarify the connection between the 
formation of radio halos and relics. 
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Notes 

1. This X-ray substructure in A2256 could have been discovered ten years before ROSAT 
since the Einstein image reveals the presence of the subcluster albeit at a lower level of 
significance (Buote 1992; Davis & Mushotzky 1993). 

2. In Buote & Tsai (1995, 1996) Hq = 80 km s"! Mpc"! and go = was assumed. 
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